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> : ABSTRACT 

a^ 

^ I We investigate the deeply virtual Compton scattering (DVCS) in the color dipole approach, 

O . implementing the dipole cross section through the saturation model, which interpolates suc- 

Q ■ cessfully between soft and hard regimes. The imaginary and real part of the DVCS amplitude 
are obtained and the results are compared to the available data. 



1 Introduction 

The investigation of hard exclusive reactions in the Bjorken limit is a reliable acces of obtaining 
information into details on the structure of the nucleon which could not be obtained considering 
inclusive deep inelastic processes. The probe provided by the photon works as a clean tool 
reactions in order to extract reliable knowledge on the substructure of strongly interacting 
particles complementar to exclusive vector meson production [1]. The recent data from the 
DESY ep collider HERA on exclusive diffractive virtual Compton process [2, 3] (DVCS) at large 
becomes an important source to study the partons, in particular gluons, inside the proton 
for nonforward kinematics and its relation with the forward one. These features are encoded 
in the formalism of the so-called skewed parton distributions, where dynamical correlations 
between partons with different momenta are taken into account. A considerable interest of 
the DVCS process lies in the particular access it gives to these generalized parton distributions 
(GDP) through the interference term with the Bethe-Heitler process. The high energy situation 
at HERA gives the important opportunity to constrain the gluon contribution to GPDs and to 
study the evolution in virtuality of the quark and gluon distributions. On the other hand. 
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Figure 1: The two leading DVCS QCD diagrams. The LO diagram in the QCD picture on the 
left graph and the NLO diagram on the right. 

recently the color dipole formalism has provided a simultaneous description of photon induced 
process. The inclusive deep inelastic reaction and the photon diffractive dissociation has been 
successfully described and the study of other exclusive processes such as DVCS is an important 
test of the color dipole approach. 

The color dipole model of diffraction is a simple unified picture of diffractive processes 
which provides a large phenomenology on the Deep Inelastic Scattering (DIS) regime. In the 
proton rest frame, the DVCS process can be seen as a succession in time of three factorisable 
subprocesses: i) the photon fluctuates in a quark- antiquark pair, ii) this coulour dipole interacts 
with the proton target, iii) the quark pair annihilates in a real photon. Using as kinematic 
variables the 7*p c.m.s. energy squared s = W'^ = {p + q)^, where p and q arc the proton and 
the photon momenta respectively, the photon virtuality squared = —q^ and the Bjorken 
scale X — Q'^/{W'^ + Q^), the corresponding amplitude at zero momentum transfer reads as, 

|*T(^,r, Q')\' 
|*L(^,r, Q')\' 

where ^'t,l are the photon wave function for transverse (T) and longitudinal (L) photons, 
respectively. They are well known from usual QED coupling, where r defines the relative 
transverse separation of the pair (dipole) and z, (l—z), are the longitudinal momentum fractions 
of the quark (antiquark). The auxiliary variable = z{l — z)Q'^ + rnj depends on the quark 
mass, rrif. The Kq^i are the McDonald functions and summation is performed over the quark 
flavors. 

The dipole cross section cr^jp describing the dipole-proton interaction is substantially af- 
fected by non-perturbative content. There are several phenomenological implementations for 
this quantity. The main feature in a successful model is to be able to match the soft (low 
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Q^) and hard (large Q^) regimes in a unified way. This has been achieved in models using 
Regge approach [4], QCD color transparency phenomena [5] and/or hybrid ones mixing non- 
perturbative functional approach and Regge phenomenology [6]. Phenomenological analysis 
including a parton saturation phenomenon at high energies has been considered, for instance 
in [7, 8]. 



In the present work we follow the quite successful saturation model [9], which interpolates 
between the small and large dipole configurations, providing color transparency behavior, (T^ip ~ 
7"^, as r ^ and constant behavior, adip ~ (Jq, at large dipolcs. The parameters of the model 
have been obtained from an adjustment to small x HERA data. Its free-parameter application 
to diffractive DIS has been also quite successful [10]. The parametrisation for the dipole cross 
section takes the eikonal-like form. 
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where the saturation scale Ql defines the onset of the saturation phenomenon, which de- 
pends on energy. The parameters were obtained from a fit to the HERA data producing 
do = 23.03(29.12) mb, A = 0.288(0.277) and xq = 3.04 • 10^^(0.41 ■ 10"^) for a 3-fiavor (4- 
flavor) analysis [9]. An additional parameter is the effective light quark mass, rrif = 0.14 GeV. 
It is worth to mention that the quark mass plays the role of a regulator for the photoproduction 
{Q"^ — 0) cross section. The charm quark mass is considered to be rUc — 1.5 GeV. A smooth 
transition to the photoproduction limit is obtained with a modification of the Bjorken variable 
as. 
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Recently, a new implementation of the model including QCD evolution [11] (labeled in the 
following BGK model) has appeared. It will be considered in the calculation to be presented 
in the last section. Now, the dipole cross section depends on the gluon distribution in the 
following way. 
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where /^^ = C/r"^ + /i^. The parameters (C and /^.q) are determined from a fit to DIS data. 
The flavor number is taken to be equal to 3. The overall normalization = 23.03 mb is 
kept flxed (labeled flt 1 in Ref. [11] )■ The function (^(a;,//^) is evolved with the leading 
order DGLAP evolution equation for the gluon density with initial scale Ql = 1 GeV^. The 
improvement preserves the main features of the low-Q^ and transition regions, while providing 
QCD evolution in the large-Q^ domain. 
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In this paper we test the vahdity of the saturation model in the particular case of the 
DVCS: 7*p — > 7p, where the final photon is real. The DVCS consists in the simplest exclusive 
diffractive process as it does not imply uncertainties due to the composite meson wave function 
in the final state like in the vector meson production, 7*p — > Vp. On the other hand, as 
wc will sec. the restriction to the transverse part of the photon wave function (at the leading 
twist), due to the real photon, enhances the contribution of larger dipole configurations and 
therefore the sensitivity, w.r.t the inclusive DIS case, to the soft content of the proton and the 
transition between the hard and soft regimes. This is why we believe that the DVCS process 
provides a particularly relevant test, in a parameter free way, of the saturation model. The 
recent cross section measurements from HI [2] and ZEUS [3] at HERA allow us now to confront 
the saturation model prediction. 

The paper is organized as follows. In the next section we address the main formulae of the 
DVCS process in the color dipole picture applying the dipole cross section given by the satura- 
tion model in order to investigate and to calculate the imaginary and real part of the referred 
amplitude. In the last section we compare the results obtained with the DVCS measurement 
of the HI Collaboration, taking into account the three(four)-flavor analysis model. The re- 
sults obtained using the recent QCD improvement to the original model are also presented. 
There, we discuss the different characteristics coming from the distinct analysis and present a 
qualitative study of the DVCS cross section using the saturation model. 

2 Deeply virtual Compton scattering in the color dipole 
picture 

In this section we establish the deeply virtual Compton scattering amplitude and the profile 
function expressions. 

Applying the color dipole approach to the DVCS process, the imaginary part of the ampli- 
tude at zero momentum transfer reads as, 

ImAis,Q'^,t^O)^ £ dz J d^r^T{z,r,Q^)adipix,r'^), (9) 
^t{z, r, = v[/^(^, r, Ql = Q^) ^t{z, r, Qj = 0) , (10) 

= ^ E 4 {[^' + (1 - ^f] ^1 ^2 Me^ r) + mj Ms, r) K,{e^ r)} , (11) 
47r f 

where 2 = z{l — z) QI 2 + f^}- 

The relative contributions from dipoles of different sizes can be analyzed with the weight 
(profile) function, 

W{r, Q2) = [' dzr ^t(z, r, Q') a,,^ {x, r^) . (12) 

J 
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Figure 2: (a) The weight function W{r,Q'^) as a function of the dipole size r at fixed c.m.s 
energy W = 75 GeV for two virtualities: = 1 GeV^ and = 10 GeV^. The solid lines 
correspond to a 3-flavor analysis whereas the long-dashed ones are for the 4^-flavor case, (b) 
Comparison between the profile function for DVCS and inclusive DIS processes at = 10 
GeV^, using the same notation. 

In Fig. (2-a), the weight function Ty(r,Q^) is shown as a function of the dipole size r at 
fixed c.m.s energy W = 75 GeV and at two virtuahties. The dipole cross section is taken 
from Eq. (4). We analyze the r-profile for a low virtuality Q"^ = 1 GeV^ and for a large one, 

= 10 GeV^. Our conclusions are similar to those in Ref. [12], i.e. as the virtualities increase 
the contribution of large dipole configurations diminishes in a sizable way. The calculation was 
performed using both a 3-flavor and a 4-fiavor analysis. The inclusion of the charm content gives 
a slightly lower normalization for the profile and by consequence for the total cross section. It is 
worth to mention that the situation is different in the inclusive DIS case, where the transverse 
wave function is symmetric on ei^2 [for instance, see Eqs. (2) and (11)] and the results including 
charm slightly differ from those using only light quarks. This is shown in Fig. (2-b) where one 
compares the DVCS and inclusive DIS profiles at = 10 GeV^ using light quarks and adding 
charm contribution. The impact of charm is smaller in the inclusive DIS case than in the DVCS 
process. Furthermore, it is verified that the DVCS profile selects larger dipole sizes in contrast 
to the inclusive DIS profile, even in the analysis including charm. Therefore, DVCS is more 
sensitive to the non-perturbative (soft) content of the scattering process. 

In order to compute the real part of the DVCS amplitude in the approach above, we have 
considered the dispersion relations as done in of Ref. [12]. In Fig. (3-a) both the imaginary 
and real part of the amplitude using the saturation model are shown as a function of c.m.s. 
energy W at = 4.5 GeV^. The ratio of the real to imaginary part, the p parameter, is 
presented in the Fig. (3-b). The analysis was done considering the light quarks (the solid lines) 
and also including the charm content (the long-dashed lines). The latter gives a lower overall 
normalization for the amplitudes, as referred before. The procedure in order to obtain the real 
part of amplitude and the p parameter was taken from Ref. [13], where a two-power adjustment 
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Figure 3: (a) The imaginary and real part of the DVCS amplitude as a function of the c.m.s. 
energy W at Q"^ = 4.5 GeV^. (b) The ratio real to imaginary part, the parameter p as a 
function of W at Q'^ = 4.5 GeV^. The solid lines correspond to the 3-flavor analysis and the 
long-dashed ones to the ^-flavor analysis. 

to the imaginary part of the amphtude is made, 

Tm^(.,t = 0)=A, (^) +A,[^] , (13) 



and the real part and p are given by, 

^ t = 0) = A, ' tan (^) + A, (^) ' tau (^^) , (14) 



0/ \^/ V'^'^o 

t2 /tj/2\A2-Ai 



tan(7rAi/2) + {W^ /W^Y''-^^ tan(7rA2/2) 

1 + A2/A1 {w^/wiY^-^^ 



(15) 



The parameters obtained from a 3-flavor (4- flavor) analysis are presented in Table (1). The 
effect in the DVCS total cross section is of order 12 % in the overall normalisation aX W = To 
GeV (p ~ 0.35). From Fig. (3-b), the p parameter depends weakly on W ioi W > 100 GeV, 
tending asymptotically to p ~ 0.3. 

Finaly, we calculate the ratio between the imaginary parts of the forward t = amplitudes 
for DIS and DVCS, R = Tm DIS/ Xm DVCS. This quantity is quite interesting as it shows the 
importance of the skewing effect. In Fig. (4) are shown the results obtained from the 3-flavor 
(solid line) and 4-flavor (dashed line) analysis of the saturation model at fixed W = 7b GeV as 
a function of the virtuality . The obtained values are slightly above those from an aligned 
jet model analysis in Ref. [14] and below those from the dipole analysis in Ref. [12]. 
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Parameter 


(3-flavor) 


(4- flavor) 


Ai 


-0.056 


-0.0336373 


A2 


0.0368 


0.0285722 


Ai 


-0.060 


-0.202454 


A2 


0.1059 


0.152216 


wi 


0.3746 


8.28475 



Table 1: Parameters of tlie two-power adjustment to tlie imaginary part of tlie amplitude, witli 
a 3-flavor (4-flavor) analysis. 
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Figure 4: The ratio of the imaginary parts of the DIS and DVCS amplitudes as a function ofQ"^ 
at fixed W = 75 GeV. The solid lines correspond to the 3-flavor analysis and the long-dashed 
ones to the ^-flavor analysis. 

3 Results and conclusions 



We are able now to compare our results with the HI Collaboration DVCS measurement at the 
photon level [2]. The final expression for the DVCS cross section is written as, 



da 

It 



1 da 
B H 



i=0 



t=o 



16 vr 



(16) 
(17) 



where B is the t slope parameter and comes from a simple exponential parametrisation and the 
p parameter has been determined in the previous section. The B value has never been measured 
for DVCS. The HI publication compares its cross section measurement to several predictions for 
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Figure 5: (a) Saturation model: a 3-flavor analysis. The photon level DVCS cross section as a 
function of the c.m.s. energy WatQ^ = 4.5 Ge\^ considering the band for the slope 5 < B < 9 
GeV'^ . (b) The DVCS cross section as a function of the virtuality at fixed W = 75 GeV, 
considering the same notation. 

5 < B < 9 GeV^^ [2]. As a first calculation, we present the results considering the saturation 
model for the 3-flavor analysis. The plots present an upper (lower) value corresponding to the 
slope values B = 5 GeV~^ [B = 9 GeV~^). In Fig. (5-a) the DVCS cross section is shown 
as a function of W at fixed = 4.5 GeV^ and in Fig. (5-b) as a function of at fixed 
W = 75 GeV. We have verified that the inclusion of charm produces a normalization slightly 
lower than using only the light quarks where one considers the same B value. This feature is 
more pronounced in the cross section as a function of energy, whereas is less significant in the 
behaviour. Concerning the choice of the quark mass values it should be noted that, taking 
m/ = increases the cross sections for the 3-flavor (4-fiavor) analysis by a factor 1.15 (1.17) at 
W=20 GeV down to 1.13 (1.16) at W=120 GeV at the average = 4.5 GeV^. This factor 
increases slightly for lower values. Keeping m/ fixed enhances the overall normalization for 
lower rUc, as usual. 

As a second calculation, we present the results considering the saturation model for the 
4-fiavor analysis. The plots present an upper (lower) value corresponding to the slope values 
B = 5 GeV^^ {B = 9 GeV~^). In the Fig. (6-a) the DVCS cross section is shown as a function 
of W at fixed = 4.5 GeV^ and in Fig. (6-b) as a function of at fixed W = 75 GeV. The 
conclusions are similar as in the 3-fiavor case. 

In what follows, due to the absence of a consistent constraint on the normalization, we will 
use the B value which fits the best with the HI measurement for each prediction. In Fig. (7) 
the following calculations are presented: (i) the saturation model for a 3-fiavor analysis (solid 
line); (ii) the saturation model for a 4-fiavor analysis (dashed line) ; (iii) the QCD improved 
saturation model (dot-dashed line), given by Eq. (8). The corresponding slope values are shown 
in the plot. We have not included the real part of the amplitude for the BGK case. The present 
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Figure 6: (a) Saturation model: a 4- flavor analysis. The photon level DVCS cross section as a 
function of the c.m.s. energy WatQ^ = 4.5 GeV^ considering the band for the slope 5 < B < 9 
GeV~^ . (b) The DVCS cross section as a function of the virtuality at fixed W = 75 GeV, 
using the same notation. 

result shows how important it is to measure the B slope in order to discriminate among the 
available approaches. It should be stressed that the result from the BGK model deserves more 
studies, since a 4-fiavor analysis was not done yet and skewed effects needed in DVCS are not 
available in this model. 

It is timely to perform a qualitative analysis in order to obtain a clear physical interpretation 
using the saturation model, as done for the DIS case in [9]. In the following we should consider 
for simplicity the particular case mj = 0. The imaginary part of the DVCS amplitude is 
given by the integrations over the dipole size r and longitudinal momentum fraction z in Eq. 
(9). One can use the small approximation argument for the Bessel function Ki{x), such that 
Ki{x) = 1/x for X and is exponentially suppressed for x ^ 1. Then, we can write: 

XmA{s,t = 0) adipix, r^) + { ] ^^dip {x, r^) (18) 

JO J4/Q2 \Q J 

where the first term is domainated by symmetric contributions, i.e. the quark and the antiquark 
constituting the dipole share almost the same momentum fraction z and the distribution on 
that variable is uniform around the mean value z = 1/2. The second term is domainated by 
aligned jet or asymmetric configurations, where z or {1 — z) takes a vanishing value and the z 
integration is restricted to z^l / {Q"^ r"^) . 

In order to proceed in the qualitative analysis, one can use an approximate form for the 
dipole cross section from the saturation model; adip = aor'^Ql/A for < 4/(5^ and adip = o"o 
for > 4/(5^. The saturation scale is given by Ql{x) = 1/Rq — x~^. Therefore, if the typical 
dipole size, l/Q, is much smaller than the mean distance between partons -Ro(^); the expression 
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Figure 7: (a) The photon level DVCS cross section as a function of virtuality at fixed W = 75 
GeV. The plots correspond to the 3-fiavor and ^-flo-vor analysis of the saturation model (solid 
and dashed lines) and the BGK model (dot-dashed line). The corresponding B slopes are also 
presented. 



for the amplitude becomes 
XmA{s,t = Q) ~ 



dj^ f apr^ Qg(x) 
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where the first term comes from symmetric configurations and the remaining ones from the 
asymmetric configurations. The final result in this case reads as, 



ImA{s,t = 0) 





1 + log 


( Q' 










1^4 A 









(20) 
(21) 



where the large logarithm in Eq. (20) comes from the intermediate region 2/Q < r < 2Rq{x) 
and the remaining comes from the regions r < 2/Q and r > 2Rq{x). The last expression gives 
the DVCS cross section at large and small x. 

The other interesting case is when the typical dipole size, 1/Q, is larger than the mean 
distance between partons, Ro{x). The contributions to the amplitude now become. 



ImA{s,t = 0) ~ / 
Jo 



dr^ / aor^gg(x) ' 

^2 \ A 



oo c?^2 

4/Q2 



g2^2 



(To, (22) 
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where the two first terms correspond to symmetric configurations and the last one to the ahgned 
jet configuration. Performing the integrals over the dipole size, the amplitude can be written 
as: 



where now the logarithmic contribution comes from the intermediate region 2Rq{x) < r < 2/Q 
and the remaining from regions r < 2Rq{x) and r > 2/Q. The last expression above gives the 
behavior of the DVCS cross section at low Q^. The qualitative results above can be contrasted 
with the DIS case, where the 7*p total cross section grows as ~ W'^^/Q'^ at high Q^, whereas 
at low it behaves as ~ In (W"^). It is worth mentioning the scaling of the amplitude in the 
ratio of the variables r = Ql{x)/Q'^, which is called the geometric scaling in the DIS case. 

In conclusion, our results are comparable with the color dipole analysis in [12], obtaining 
a good agreement with DVCS cross sections on the W and dependences. We have tested 
the saturation model in the exclusive DVCS process, showing its consistency. The utilized 
dipole cross section interpolates successfully the soft and hard regions, producing a smooth 
transition between them. In particular, the sensibility to the inclusion of charm content seems 
to be more pronounced compared with the inclusive DIS case. This can be investigated in more 
details when measurements with higher virtuality will become available. A simple qualitative 
picture for the DVCS process can be drawn, presenting a clear scaling property. The QCD 
improvement provides good agreement with data requiring a higher B slope in contrast with 
the original model. The present analysis also shows the importance of a measurement of the B 
slope for the DVCS process. 
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